Example of WISP.data functions
Alessandro Oggioni
Source:vignettes/exampleOfFunctions.Rmd
exampleOfFunctions.RmdProvide the account of your station
userid <- Sys.getenv("WISP_IREA_USERID")
pwd <- Sys.getenv("WISP_IREA_PSW")Get data from the WISP station server for a multiple dates
If you want to download the data for a single date, consider using
the function wisp_get_reflectance_data() instead.
reflect_data <- WISP.data::wisp_get_reflectance_multi_data(
time_from = "2024-04-08T09:00",
time_to = "2024-04-10T14:00",
station = "WISPstation012",
userid = userid,
pwd = pwd
)Quality Check (QC) of the reflectance data
reflect_data_qc <- WISP.data::wisp_qc_reflectance_data(
data = reflect_data,
maxPeak = 0.05
)##
## ----
## 9 spectral signatures were removed during QC:
## - 5 spectral signatures were removed thanks to QC1
## - 4 spectral signatures were removed thanks to QC2
##
## The spectral signature of 2024-04-08 10:45:05 has been removed thanks to: QC2
## The spectral signature of 2024-04-08 11:00:05 has been removed thanks to: QC2
## The spectral signature of 2024-04-08 11:15:05 has been removed thanks to: QC2
## The spectral signature of 2024-04-08 11:30:05 has been removed thanks to: QC2
## The spectral signature of 2024-04-08 13:00:04 has been removed thanks to: QC1
## The spectral signature of 2024-04-08 13:15:05 has been removed thanks to: QC1
## The spectral signature of 2024-04-08 13:30:05 has been removed thanks to: QC1
## The spectral signature of 2024-04-08 13:45:04 has been removed thanks to: QC1
## The spectral signature of 2024-04-10 13:45:05 has been removed thanks to: QC1
##
## QC1 remove spectral signatures with negative values below 845 nm
## QC2 remove spectral signatures with outliers in the NIR (840 nm > 700 nm)
## ----
Sunglint Removal (SR) of the reflectance data
In this function if the save_output parameter is set to
TRUE, the function will save the reflectance data with the
sunglint removed as a .csv file in the working directory.
reflect_data_sr <- WISP.data::wisp_sr_reflectance_data(
qc_data = reflect_data_qc,
save_output = TRUE
)Plot the reflectance data
WISP.data::wisp_plot_reflectance_data(
data = reflect_data_sr,
legend_TSM = TRUE,
legend_Chla = TRUE,
legend_Kd = FALSE,
legend_cpc = FALSE
)